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TRANSPORT COEFFICIENTS IN ROTATING WEAKLY COMPRESSIBLE 

TURBULENCE * 


ROBERT RUBINSTEIN*, YE ZHOU*, AND GORDON ERLEBACHER 5 

Abstract. Analytical studies of compressible turbulence have found that compressible velocity fluctua- 
tions create both effective fluid transport properties and an effective equation of state. This paper investigates 
the effects of rotation on compressible turbulence. It is shown that rotation modifies the transport proper- 
ties of compressible turbulence by replacing the turbulence time scale by a rotational time scale, much as 
rotation modifies the transport properties of incompressible turbulence. But thermal equilibrium properties 
are modified in a more complex manner. Two regimes are possible: one dominated by incompressible fluc- 
tuations, in which the sound speed is modified as it is in non-rotating turbulence, and a rotation dominated 
regime in which the sound speed enhancement is rotation dependent. The dimensionless parameter which 
discriminates between regimes is identified. In general, rotation is found to suppress the effects of compress- 
ibility. A novel feature of the present analysis is the use of a non- Kolmogorov steady state as the reference 
state of turbulence, introduction of such steady states expands the power and utility of analytical turbulence 
closures to a wider range of problems. 

Key words, compressible turbulence, rotating turbulence, closure models 

Subject classification. Fluid Mechanics 

1. Introduction. Rotating weakly compressible turbulence occurs in astrophysical flows and in engi- 
neering flows like swirling jets. The development of transport models and sub-grid scale models for such 
flows is difficult because the coupling between compressibility and rotation introduces at least two natural 
dimensionless parameters: the turbulent Mach number Mt = K l ^/c and the inverse turbulent Rossby num- 
ber ( = f IK/e. It is not easy to assess the combined effects of these parameters heuristically. Recourse to 
DNS will require an unrealistically large number of simulations to cover a useful range of parameters. LES 
will not be practical without sub-grid scale models which are sensitive to both rotation and compressibility: 
the development of such models is a goal of this paper. 

The present work addresses these problems using the rational statistical turbulence closure provided by 
Yoshizawa’s (1986,1995) two-scale direct interaction approximation (TSDIA), which is a natural and analyt- 
ically tractable extension of Kraichnan’s (1959) direct interaction approximation (DIA) to inhomogeneous 
turbulence. It is based on previous work on the transport theory of weakly compressible turbulence (Rubin- 
stein and Erlebacher, 1997, RE in what follows) and on an analytical theory of rotating turbulence initiated 
by Zhou (1995). 

The problem of weakly compressible turbulence, in which the equations for compressible quantities are 
linearized about a basic state of incompressible turbulence, was considered from the standpoint of TSDIA 
in RE. TSDIA leads to a description of compressibility effects on the mean flow equations by a set of trans- 
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port coefficients: this description generalizes the eddy viscosity description of incompressible turbulence. 
Following Yoshizawa, the transport coefficients are expressed in terms of the fundamental statistical descrip- 
tors of weakly compressible turbulence: the response and correlation functions of incompressible turbulence 
(Kraichnan, 1959) and the additional statistical descriptors of the compressible part of the motion which 
will be defined later in this paper. By assuming that the basic state of incompressible turbulence is in a 
Kolmogorov steady state, it is possible to derive both formulas for the transport coefficients in terms of single 
point descriptors of turbulence and, following Yakhot and Orszag (1986), sub- grid scale models for LES of 
weakly compressible turbulence. 

The analysis of rotation effects will be based on Zhou’s (1995) scaling theory of rotating turbulence, 
which demonstrates the existence of a steady state with constant intertial range energy flux in rotating 
turbulence, but with non- Kolmogorov scaling exponents. In the limit of strong rotation, this theory is closely 
related to the theory of weak turbulence (Zakharov et al } 1992) in which nonlinear temporal decorrelation of 
Fourier modes is dominated by linear dispersive decorrelation. The extension of the theory of RE to rotating 
turbulence will be accomplished by perturbing about this state of rotating incompressible turbulence instead 
of about the Kolmogorov state. 

Like previous investigations of compressible turbulence (Staroselsky et ai, 1990; Chandrasekhar, 1951), 
RE found two kinds of compressibility effects: modified transport coefficients and modified equilibrium 
properties. Turbulence creates diffusivities for entropy and pressure and enhances the molecular shear and 
bulk viscosities. But turbulence also modifies the specific heat ratio and creates an effective turbulent 
pressure. These combined effects lead to an enhanced propagation speed for sound waves in the mean 
flow; they can be considered as a modified equation of state for compressible turbulence. Rotation alters 
the transport properties, as it does in incompressible turbulence, by replacing the turbulent time scale by 
a rotational time scale. It will be shown that rotating compressible turbulence can exhibit two different 
effective equations of state, depending on the size of the parameter In general, rotation tends to 

suppress or counteract the effects of compressibility. 


2. Review of the TSDIA Theory of Weakly Compressible Turbulence. The basic result of RE 
is the following model for weakly compressible turbulence: 


(1) 

( 2 ) 
( 3 ) 


— + U • vs = V • (u ss VS) 

™ + u • vu + (i - 4 )VP = V • {-\ ki + *[vu + (vum 


dP 

dt 


+ V(i/“ U V • U) 

+ U • VP + (7 - 7)PV • U = V • (t/^VP) 


The turbulence dependent properties are defined by the integrals 
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The quantities ia i/ 3S , i/ pp and t/ uu are effective transport properties, whereas the effective turbulent pressure 
n and effective mean density R define the modified equation of state of compressible turbulence. In Eqs. 
(4)-(10), G and Q are the response and correlation functions, the fundamental statistical descriptors of 
isotropic turbulence in DIA (Kraichnan, 1959) and Q v is the twotime pressure correlation function. The 
compressible response functions G * , G** , G ™ , G ”* , G™ are defined in terms of the normalized Helmholtz 

decomposition 


iii(k, i) = ^t(k, t) + ih{k ^>(k, t) 

and satisfy the DIA response equations, written in terms of time Fourier transforms, 

(12) iwG** (k, w) + t] ww G ww (k, u>) = 1 

and 
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Following RE, the response equations will be solved in a Markovianized limit in which the damping functions 
rj in Eqs. (12), (13) are defined by 

T7 w “(k) = --Py(k) [ dr [ dpdqrJ mn (k,p,q)G ,m "(q,T) x 

2 do dk=p+q 

— P? k)Q mr (p, t) 

- - Pij (k) / dr f dpdqT j mn (k, p, q)G** (q, r) x 

2 Jo Jk=p+q 
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( 18 ) 


*/*<& = - dr dpd q r^ l (k, p,q)G rPP (q, r) x 

»/o »/k=p-fq 

r*(q, -p,k)<? mn (p,r) 

where the quantities T, which characterize the interactions between the fields, are given by 

P> q) ~ i{Qm&in "I" Pn^im) 
r n(k,p,q) = ik n 

These equations are solved approximately in RE by a series expansion in turbulent Mach number. Weak 
compressibility assumes the existence of a basic state of incompressible turbulence. If this turbulence is in 
a Kolmogorov steady state, in which isotropy implies that 

( 19 ) G«(k,r)=G(fc,r)Py(k) 

< 2 °) Qij(Kr)=Q(k,r)P ij ( k) 

the approximation leads, at lowest nontrivial order, to the conclusion that G ww is the response function of 
a passive vector field, and to the expressions 


(21) 

G* p (k,T) = -^exp(-r 7 ^r/2)sin(5r)//(r) 
. £>2 

G p f(k, r) = p 0 — exp(-r7^T/2)sin(5r)ff(r) 
G pp (k,r) = exp(— t7^’t/2)cos(<St)//'(t) 

+ 7^ exp (-r^ r /2) sm(Sr)H(r) 


G^(k,T) = exp(-^r/2)cos(5r)jy(r) 

(22) 

rj^(k) 

2< 5 ; exp(-^r/2)sin(5r)J/(r) 

where 


(23) 

5 = { c 2 k 2 - ^^(A:) 2 } 1 / 2 


In non-rotating turbulence, the damping function rf* is determined by the Reynolds analogy 


< 24 ) V*+(k) = a^ik) 

where a# is an inverse Prandtl number, and the incompressible spectral damping function rj(k) takes the 
Kolmogorov scaling form 


( 25 ) r,(k) = C D e ^k 2 ' 3 

Analogs of these formulas for rotating turbulence will be developed subsequently. Once such formulas 
are known, the integrals for transport coefficients Eqs. (4)-(10) can be evaluated explicitly. For turbulence 
transport models, the integration will extend over the entire inertial range k > k 0 where k 0 is the inverse 
integral scale. Subgrid models are derived instead (Yakhot and Orszag, 1986) by integrating only over 
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Table 1 

Scaling results of RE 



Single-point model 

General subgrid Model 

Smagorinsky model 

v vv 

M?K 2 /£ 

eA 2 /c 2 

(5 < ) 3 A 4 /c 2 

U uu 

K 2 /e 

£ l/3 A 4/3 

(S<) A 2 

n 

KM 2 

£ 4/3 A 4/3 /c 2 

(5<) 4 A 4 /c 2 

l/R 

M? 

£ 2/3 A 2/3 /c 2 

(5<) 2 A 2 /c 2 


wavevectors which satisfy the condition k > 2n/A y where A is the filter size; as usual, this derivation assumes 
that these scales are in the inertial range. The quantity e can be closed by Smagorinsky’s hypothesis of local 
energy equilibrium. An interesting alternative, related to postulating the sweeping hypothesis for turbulent 
time correlations, has been suggested recently by Yoshizawa et al (1996). 

The results of RE are summarized in Table 1. The effective properties are given without the appropriate 
constants, all of which could in principle be evaluated theoretically. The column labelled “general subgrid 
moder contains the result of integrating the integral expression over the subgrid scales; the corresponding 
“Smagorinsky model” is obtained by equating the dissipation rate to the resolved production. Throughout 
the table, S K is defined by (S < ) 2 = S <- , where is the resolved fluctuating strain rate. 

3. Review of the Statistical Theory of Rotating Turbulence. Kolmogorov scaling applies to a 
homogeneous turbulent steady state with constant energy flux. External effects like rotation will modify 
this steady state. Rotation is an especially simple external effect, since energy remains an inviscid invariant 
under rotation, and a steady state with constant energy flux remains possible. Kolmogorov scaling no longer 
applies to this steady state and the existence of a distinguished time scale, namely the inverse rotation rate, 
precludes the deduction of the applicable scaling law by dimensional analysis alone. To deduce the scaling, 
we will appeal to closure in the form of the direct interaction approximation. 

For rotating turbulence, the DIA equations of motion take the form 


(26) 


G i:j (k, t, s) + 2P ip {k)n pq G qi (k, t , s) 

+ J dr T) ip (k,t,r)G pj (k,r, s) = 0 
Qij(k,t,s) + 2P lp {k)Cl pg Q qj (k,t, s) 

+ J dr r? ip (k, t, r)Q pj (k, r, 5 ) 

= / dr G ip (k,t,r)F pj (k, s,r) 

Jo 

where the eddy damping tj and forcing F are defined by 

T) ir (k,t,s)= dpdt{Pi Tnn (k)P flra (p) x 

J k=p+q 


(27) 
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(28) 


Cm/i (P ,t, s)Q ns (q, t, s) 

F l} (k, t, s) = f dpdqP imn (k)P,> s (k) x 

J k=p+q 

(29) Qns(pj ^ 

In Eqs. (28) and (29), 


Ptmn(k) = k m P in {k) + k n Pim{ k) 
P%m (h) = ^im k k{ k m 


and £l p<z is the antisymmetric rotation matrix. 

The solution of these equations in complete generality is not known. A natural perturbation theory treats 
the rotation terms as small, and perturbs about an isotropic turbulent state. This approach is adopted by 
Shimomura and Yoshizawa (1986), who derive a TSDIA theory in which inhomogeneity and rotation are 
both described by small parameters. 

A complementary limit is also of interest. Namely, in the response equation, balance the time derivative 
by the rotation term, and treat the eddy damping as small. This linear theory of the response equation 
treats strongly rotating turbulence as a case of weak turbulence (Zakharov et al , 1992) in which nonlinear 
decorrelation of Fourier modes is dominated by linear dispersive decorrelation (Waleffe, 1993). The result is 
conveniently expressed in terms of the Craya-Herring basis 

e (1) (k) = kxO/ | kxfl \ 
e< 2 >(k) = kx(kxft)/ | kx(kxO) | 

or the equivalent helical mode basis used, for example, by Waleffe (1993) and Cambon and Jacquin (1989), 
and the corresponding tensors 

f?. 

Stj 

f 2 

(30) 

Note that (Jj = Pij (k) . 

In this approximation, the leading order solution of Eq. (26) is 

(31) Gij( k, r) = {cos(0r)Pij(k) + sin(0T)£° (k)}tf (r) 

In Eq. (31), 0(k) is defined by 


- P (i)pp) _ 

= e U) e (2) +e (D e ( 2 ) 
_ e (l) e (l) _ e (2)J2) 

— e i e- 

= ej 1 >e< 1) +e< 2) ef 


0(k) = k«/fc 

r = t — s is the time difference, and H is the unit step function. 

Let us adopt the fluctuation- dissipation hypothesis relating the two-time correlation function to the 
response function and single-time correlation function 

(32) Qij( k,r) = Gj m (k, r)Q mj (k) + G jm (k, -r)Q mi (k) 
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Conditions under which this approximation is reasonable are discussed by Woodruff (1992, 1994). Substitut 
ing Eqs. (31) and (32) in Eq. (27) suggests that the single-time correlation function should take the general 
form containing all of the £ tensors of Eq. (30) 

(33) Qa( k) = £ Q p ( k)^-(k) 

0<p<3 

which is equivalent to the form of the correlation function noted by Cambon and Jacquin (1989). But 
substituting Eq. (33) in Eq. (32), we find instead that the symmetry 

Qij{ k,r) = Qji( k, -r) 

leads at lowest order in the weak turbulence approximation to the simpler expression 

(34) Qij (k, r) = sin(©T)Q(k)£° (k) + cos(0r )Q( k)^ (k) 

This simplification is a consequence of the approximate treatment of time correlations, refinement of the 
approximation will recover the full expression Eq. (33). 

The DIA inertial range energy balance (Kraichnan, 1971), which states that a steady state with constant 

energy flux exists, is 

rOC 

e = [1^ — I ]f > t mn (k) j dx 2P^ rs (p)G m/x (p, r)Q ns (q, (k, r) 

Jo 

(35) - P jrs (k)Gij(k , r)Q ns (p, r)Q mr { q, r) 

where the integration operators in Eq. (35) are defined by 

J+(fc 0 )= f dk I <5(k-p-q)dpdq 

J k>ko Jp,q<k o 

I~{k 0 )~ j dk j <5(k - p - q)dpdq 

Jk<ko J o 

Substituting Eqs. (31) and (34) in Eq. (35) shows that the time integrals in Eq. (35) will contain the factors 

(36) T(k, p, q) = fl _1 <5(±p • fl/pfl ± q • Cl/qSl ± k • il/kQ) 
where 

n = (n pq n pq ) 1/2 

Thus, wave-vector integrations take place over resonant triads only (Waleffe, 1993) and these integrals scale 
as SI -1 (Zhou, 1995). 

Following the derivation of the Kolmogorov steady state, we seek a solution of Eq. (35) of the power- law 
form 

(37) Q( k) = k-°- 2 f(n ■ k /Ml) 

In view of Eq. (36), the energy balance Eq. (35) requires the scaling a = 2 and the proportionality / ~ \feVt- 
Define the shell-averaged energy spectrum 


(38) 


E(k) = 
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then (Zhou, 1995) 


(39) 


E(k) = C%VdHk~ 2 


It is easily verified that the flux integral Eq. (35) converges for an infinite k~ 2 inertial range. Consequently, 
the scaling law of Eq. (39) is the scaling of a solution of the DIA equations of motion in the weak turbulence 
approximation . 

The anisotropy of rotating turbulence implies that the eddy viscosity is a tensor defined so that 

d U r 

< UiUj > = Vijra 

where TSDIA gives 


( 40 ) «^p. = \ J dk j dr G ir (k, r)Q js { k, r) + ( ij)(rs ) 

and the parentheses indicate index symmetrization. Using the general form of the single-time correlation 
function Eq. (33) and performing the time integration, 

(4!) v ijTS = Jdk i<5(: k • n/kn) ^^(k)C( k )Q A ( k ) 

The delta function dependence of the integrand of Eq. (41) implies that this integrand is nonvanishing 
only for vectors perpendicular to the rotation axis. This fact originates in the resonance condition noted 
earlier: since TSDIA treats the effect of small turbulent scales on large scales characterizing the mean flow, 
it evaluates turbulence transport properties in the distant interaction approximation (Kraichnan, 1987), in 
which interactions among modes with wavevectors k, p, q are all treated as distant interactions for which 
k/p, k/q — * 0. In this limit, only planar triads such that k fl = p fl = q fl = 0 can satisfy the resonance 
condition. This states the “two-dimensionalization” of energy transfer in strongly rotating turbulence, an 
observation which is consistent with the more comprehensive theory of Mahalov and Zhou (1996). 

The isotropic contribution to the eddy viscosity, obtained by multiplying Eq. (41) by the isotropic tensor 
fiijdrs + 5i r djs + Si 9 Sj r , is the scalar 


<«) -P?5 

For a subgrid- scale model, instead 


(43) i/(A) = Cgyfi/QA 
Closing the dissipation rate e by the resolved production 

(44) e = ^(A)5^ = !/(A)(5<) 2 

where is the resolved fluctuating strain rate, leads to the Smagorinsky model for strongly rotating turbu- 
lence 


(45) i/(A) = C2 2 {S<) 2 A 2 /Q 

Thus, the subgrid scale viscosity is suppressed by rotation. The eddy viscosity has the same A 2 dependence 
on grid size as the standard Smagorinsky model for non-rotating turbulence. 
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Unlike a fully linear theory like rapid distortion theory, the weak turbulence approximation ignores 
nonlinearity only in the response equation but retains the nonlinearity of the energy flux condition. The 
corrections due to nonlinearity of the response equation can be systematically evaluated in the weak turbu- 
lence approximation. Namely, substituting Eqs. (31) and (39) in Eq. (28) for the eddy damping factor and 
taking the distant interaction limit, 


(46) 


roo | 

k 2 / dp —\fsVtp~ 2 ~ ky/ejft 

Jk “ 


The corrections to the time scale and energy spectrum therefore have the form 


(47) 


t ~ ^{i + o(fr 3/2 )} _1 

E ~ \feQ.k~ 2 { 1 + 0((k 2 e/fl 3 ) 1/2 } 


The low rotation rate expansion of Shimomura and Yoshizawa (1986) gives the complementary expansions 
in positive powers of ft, 

T~e~ 1/3 fc“ 2/3 {l + 0(H)} 

(48) E ~ e 2 ' 3 k~ 5 ' 3 { 1 + 0(n/e'' 3 k 2 ' 3 } 

4. The Effect of Rotation on Weakly Compressible Turbulence. In the presence of rotation, 
centrifugal and Coriolis terms must be added to the momentum equations. However, the effects of rotation on 
the local thermal equilibrium properties of fluids can usually be neglected, so that the entropy and pressure 
equations remain unchanged. The modified momentum equations have the form 

(49) Ui + 2n tp u p -{l + fyn 2 Xi = ■■■ 

where terms independent of rotation have not been explicitly written. The mean contribution to the centrifu- 
gal term should be added to the mean momentum equation. The fluctuating component has both solenoidal 
and potential parts, which will appear in the equations for and <f>. 

The Coriolis term must be treated by the Helmholtz decomposition Eq. (12) followed by a second 
decomposition into solenoidal and potential parts, 


n ip u p = [P tJ (k) + P,*(k)]U JP [u* +w p + ik~ l k p 4>\ 

— Pij{ k)fi, p u£° + PipSlpjWp + ik~ 1 il ip k p (p 

(50) + Pi p ^jp u p T PipQpj w p 

The first term in Eq. (50) is the usual Coriolis term which causes the modified statistically steady state 
described in Sect. 3. The second term is the analogous Coriolis term for the compressible field The 
third term is transverse; therefore, it couples the field <f> to the Wi equation. This coupling does not occur in 
non-rotating turbulence. The fourth term shows that the incompressible field is a source for the compressible 
potential field <p. Finally, the fifth term couples the field w l to the (f) equation. The anisotropy introduced 
by rotation makes this type of coupling between a scalar and a vector possible; this coupling is obviously 
impossible in isotropic turbulence. Note the important fact that <f> is not coupled linearly to itself through 
rotation. 

In the analysis of non-rotating weakly compressible turbulence, it was not possible to treat all couplings 
exactly; likewise, the present analysis will treat all rotation- dependent couplings perturbatively, except for 
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the rotational self-couplings of the solenoidal fields w and u°°. This defines the diagonal approximation for 
rotational effects 


u°° + 2P ij {l jp u™ = ■ ■ ■ 
ibi -f- — iJc ~\~ 

(51) 0 = -2f>;n jp (u- + u; p )... 

where dots indicate centrifugal and rotation- independent terms. In this approximation, it is appropriate to 
replace the Kolmogorov steady state for the basic incompressible field u°° by the k ~ 2 steady state for strong 
rotation described in Sect. 3. 

5. Evaluation of the Transport Coefficients. 

5.1. The entropy diffusivity. Since the entropy equation is not altered by rotation, G ss remains the 
response function of a passive scalar. The usual Reynolds analogy, which postulates the proportionality of 
the time correlations of the velocity and passive scalar fields leads, following the derivation of Eq. (42) to 

(52) 

in the strong rotation limit. 

5.2. The quantities II and R . According to Eqs. (8) and (10), the evaluation of these quantities 
requires the response function G^ p and the two-time pressure correlation in rotating incompressible turbu- 
lence. The perturbation theory of RE begins by observing that G ww is the response function of a passive 
vector field. Whereas this could be asserted exactly in non-rotating turbulence, the couplings which occur 
in Eq. (51) restrict this assertion to the lowest order in the diagonal approximation for rotation effects. 

Thus, we assume that to lowest order G ww takes the same form as the response function in Eq. (31) 

(53) G%°( k, r) = {cos(0r)P ii (k) + sin(er)£° (k)}tf (r) 

As in RE, the damping functions except for all vanish to this order. We must first check that the integral 
in Eq. (15) which defines this function converges for the k~ 2 spectral scaling. As usual, this convergence 
check requires the evaluation of the geometric factors in this integral in the limits p, q — ► oo. This factor has 
the form 


(<?m^in + kmfiin)(krdjp Pj^pr)Cmr (p)Cnp (q) 

= fcmfcr£mr(P)f£ (q) ~ *mMmp(P)£,>(q) 

(54) + fc n fcrf,y(p)S£,(q) - *nM,>(p)££ p (q) 

The second and fourth terms, which power counting would suggest cause logarithmic divergence in Eq. (15), 
instead vanish by symmetry. The remaining terms are of order p~ 2 and therefore converge for large p. The 
infrared divergence when p — *■ 0 is explained by the sweeping effect on Eulerian time correlations to which 
Kraichnan (1964) first called attention. Rather than evaluate the transport coefficients in a Lagrangian 
statistical theory (Kraichnan, 1965; Kaneda, 1981) in which this divergence does not occur, we shall sim- 
ply regularize the Eulerian theory by integrating over the range p>q > k only (Kraichnan, 1963). This 
regularization leads finally to the result 

(55) »?*(k,T) = c&VS7n* 
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which can be compared to the correction term computed in Eq. (46). 

The evaluation of the two-time pressure correlation function is actually simpler in the weak turbulence 
approximation for rotating turbulence that it is for non-rotating turbulence because the resonance condition 
and even time parity of time correlation functions force the time dependence cos(0t). As in the non-rotating 
case, convergence of the single- time correlation under the quasi- normal approximation appropriate to the 
direct interaction approximation is easily established, and we conclude that 

(56) Q p (k , r) = C' J B,ncos(0r )eQk~ 3 / Airk 2 


It is now straightforward to verily the calculation 


(57) 

where 


n/7 = ~ J dk J dr exip(—V0kr) x 

[sin^/c + Q.)t + sin(V r /c - Q)r]Q p (k) 


f 


vk 3 {o) + v 2 )-vkn 2 


’ „ eQ w 9 

Arrrl ic 9 h~ Z 

Z V [0 + V 2 ] 2 k 4 + 2(0 - V 2 )Q 2 k 2 + n 4 


V = v/c 2 -6» 


and 


The integral in Eq. (57) could be evaluated in closed form. However, its limits for large and small wavenumber 
k are easily evaluated directly. Assume that the k~ 2 spectrum will be cut off at wavenumber k 0i which is 
either the inverse integral scale, or inverse filter size. If the cutoff scale /cq is large, then approximately 


(58) 


n/7 = 4n 



-eflk 2 

c 



r-2 


Setting ko to the inverse integral scale leads to the single-point model 
(59) n/7 = Cg>^- 

These results coincide with those of RE except that the energy spectrum and total kinetic energy are 
determined from the k~ 2 spectrum instead of from the Kolmogorov spectrum. However, the large scales 
contribute a different limit, namely, 


fk* UJL _ 

(60) n/7 = 4tt / —petlK- 2 = -log(fcVfco) 

J fco 

The scale k* is the scale at which the high wavenumber limit Eq. (58) applies; its value is calculated below. 
Note the logarithmic contribution. 

The approximation Eq. (58) applies to small scales which are relatively insensitive to rotation. Whereas 
the theory of weak compressibility requires the condition K << <? , the limit Eq. (58) requires the much 
stronger condition 


(61) 


K « cyfijn 
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At very high rotation rates, a range of large scales such that >> ck can exist. Such scales are described 
instead by the limit Eq. (60). The transition between spectral regions is determined by the parameter 

(62) r = - 

c 

Since Eq. (60) describes a condition in which compressibility effects have been reduced, we conclude that 
rotation tends to counteract the effects of compressibility. The single-point parameter which discriminates 
between regimes is the ratio 

is 1^1/2 

(63) _^ = ^L_ v ^ =c i/2 M< 

Cyje/U C 

Weak compressibility is defined by the condition that M t is small. When in addition i s small, rotation 

effects are weak; strong rotation effects occur only if is large. Note that since rotation reduces the 

energy transfer £, large values of the parameter C are possible in rotating compressible turbulence. 

The evaluation of the coefficient R proceeds similarly. However, it should be noted that the anisotropy 
of rotating turbulence actually makes this quantity a tensor, so that the mean equations of motion contain 
the pressure term 




1 , dP 

~ 0^17 


R t3l dX } 


instead of the isotropic pressure term of Eq. (3). The possibility that sound waves in the mean velocity 
field become dispersive as a result of the direction dependence of the effective sound speed deserves further 
investigation. In this article, the sound speed will be treated as a scalar. By comparing Eqs. (8) and (10), it 
is evident that R is found from the previous calculation by multiplying the integrand in Eq. (8) by ky/F/Tl. 
The small scale limit becomes 


(64) 

and the large scale limit becomes 


1 _ s-in eQk 0 z &o 

k~° R> 


2 L Ve^o 1 

Veto R> c 2 


(65) R Cr< 

The single-point quantities corresponding to Eqs. (64) and (65) are simply 

Af f 2 for /co — > oc 
1/C for k 0 — > 0 

Eq. (2) shows that R represents a modification of the mean density due to compressible flucations. Eqs. 
(65) and (66) show that this modified mean density actually vanishes for very large rotation rates, again 
suggesting that rotation tends to counteract compressibility effects. 



5*3. The pressure and bulk viscosities. Like R , these quantities are tensors in rotating turbulence, 
but only scalar values will be computed here. Eq. (6) shows that the pressure diffusivity v pp is the sum of 
two terms. For the first term, in the low turbulent Mach number limit, 




ik. 


dk 


if /A + 2i)S 2 


(67) 


ko (r,/ 2) 4 + 2(r?/2) 2 (5 2 + fi 2 ) + (S 2 - Q 2 )2 

roc 

/ dk \/7/Ct(c 2 - 36/4)~ 1 k~ 1 VeQk~ 2 
Jk 0 

K 2 


VeQk 2 


tho 

= C r 


fi(c 2 -3(9/4) 
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Table 2 

Main Scaling Results 



Single-point model 

General subgrid Model 

Smagorinsky model 

v vv 

M?K/Q 

£ A 2 /c? 

(5<) 3 A 4 /c 2 

V uu 

K/n 

^/TjUA 

(5<) 2 A 2 /fl 

II (large k) 

KM 2 

eflA 2 /c 2 

n(S ,< ) 3 A 4 /c 2 

11 (small k) 

e/n 

eA°/fl 

(S<) 3 A 2 /fi 

l/R (large k) 

M 2 

v/rfiA/c 2 

(S < ) 3 / 2 f} 1 / 2 A 2 /c 2 

l/R (small k) 

e/KCl 

yje/ SPA- 1 

(5<) 3 / 2 n^ 3 / 2 A° 


where 77 is the incompressible damping factor of Eq. (25) This expression applies to both large and small 
values of the cutoff scale Icq. The second term is evaluated like the entropy diffusivity as 


( 68 ) 




i* p = ^fc - 2 = C 2 


AT 2 

He 2 


In summing these contributions, care must be taken that cancellations between these terms do not change 
the limits of the transport coefficients. In the absence of such cancellations, 


(69) 


,,pp 


pp c 2 n 


This result can be compared with the non-rotating result in RE. Rotation has the same effect that it has on 
the incompressible eddy viscosity: the turbulence time scale K/e is replaced by 1/S7. 

A similar calculation leads to an explicit expression for the bulk viscosity, 


(70) yUU = C n% 

Note that in RE, this quantity was of order M t °. Again, in the limit of strong rotation, the consequence of 
rotation is the simple replacement of the incompressible turbulence time scale by 1 /fb 

The results of this section are summarized in Table 2. The format follows Table 1, except that both 
large and small scale limits are given for the equilibrium properties It and R. 

6. The Equation of State of Rotating Weakly Compressible Turbulence. Beginning with the 
first mixing length models, turbulence modeling has always postulated that turbulent fluctuations modify 
fluid transport properties by replacing thermal fluctuations with scale dependent macroscopic turbulent 
fluctuations. Analyses of compressible turbulence (Chandrasekhar, 1952, Staroselsky et a/, 1990) have noted 
analogous turbulence effects on thermal equilibrium properties: velocity fluctuations lead to a modified 
equation of state of compressible turbulence. Thus, Eqs. ( 6 ), ( 8 ) and (10) show that fluctuations modify 7 
and lead to changes of the effective pressure and density represented by the terms II and R. 
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In rotating compressible turbulence, the picture is somewhat more complex. The transport coefficients 
v are all modified like the incompressible eddy viscosity: the turbulence time scale K/e is replaced by the 
rotation time scale l/f2. But thermal equilibrium properties are modified so that two different regimes are 
possible: a regime dominated by incompressible fluctuations, in which the sound speed is modified as it is in 
non-rotating turbulence, and a rotation dominated regime in which the sound speed enhancement is propor- 
tional to \J e/Q instead. As noted earlier, this limit indicates a reduction of compressibility effects on large 
scales by rotation. Pure compressibility effects enter the evolution equations through the turbulent transport 
and thermal equilibrium properties of Eqs. (2)- (4) which are evaluated in Table 1; the further effects of 
rotation on the evolution equation appear through the rotation-induced modifications of the transport and 
equilibrium properties which are evaluated in Table 2. 

7 . Conclusions. Yoshizawa’s TSDIA formalism provides a systematic approach to the complex coupled 
field problem of rotating weakly compressible turbulence. In this problem, the appearance of two natural 
dimensionless parameters would frustrate attempts to derive appropriate models heuristically. In particular, 
the modified equation of state in rotating weakly compressible turbulence is not obvious. 

The main defect of this calculation is the assumption of a diagonal approximation for rotation effects, 
although like RE, we do not make a diagonal approximation for compressibility effects. We should attempt 
either to verify that the non-diagonal terms are small, or look for a more sophisticated approximation. 

A novel feature of the present analysis is the use of a non- Kolmogorov steady state as the reference state 
of turbulence in a TSDIA calculation. The introduction of such steady states expands the power and utility 
of TSDIA to a wider range of problems. 

Compressible swirling jets are a natural class of flows to which the present theory should apply. The 
theory is consistent with the general observation that rotation counteracts the suppression of mixing which 
is the physically most prominent feature of compressible jets. But quantitative conclusions will obviously 
require computational tests of the theory. 
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